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Abstract: We prove a comparison inequality between a system of indepen- 
dent random walkers and a system of random walkers which either interact 
by attracting each other - a process which we call here the symmetric in- 
elusion process (SIP) - or repel each other - a generalized version of the 
Q_( . well-known symmetric exclusion process. As an application, new correla- 

r-| ! tion inequalities are obtained for the SIP, as well as for some interacting 

diffusions which are used as models of heat conduction, - the so-called 
Brownian momentum process, and the Brownian energy process. These 
inequalities are counterparts of the inequalities (in the opposite direction) 
for the symmetric exclusion process, showing that the SIP is a natural 
bosonic analogue of the symmetric exclusion process, which is fermionic. 
\Q \ Finally, we consider a boundary driven version of the SIP for which we 

| prove duality and then obtain correlation inequalities. 
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1 Introduction 

In Liggett |14j . Chapter VIII, proposition 1.7, a comparison inequality between inde- 
pendent symmetric random walkers and corresponding exclusion symmetric random 
walkers is obtained. This inequality plays a crucial role in the understanding of the 
exclusion process (SEP); it makes rigorous the intuitive picture that symmetric ran- 
dom walkers interacting by exclusion are more spread out than the corresponding 
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independent walkers, as a consequence of their repulsive interaction (exclusion), or 
in more physical terms, because of the fermionic nature of the exclusion process. 
The comparison inequality is a key ingredient in the ergodic theory of the symmet- 
ric exclusion process, i.e., in the characterization of the invariant measures, and the 
measures which are in the course of time attracted to a given invariant measure. The 
comparison inequality has been generalized later on by Andjel [I], Liggett [T5] , and 
recently in the work of Borcea, Branden and Liggett [3] . 

In the search of a natural conservative particle system where the opposite inequal- 
ity holds, i.e., where the particles are less spread out than corresponding independent 
random walkers, it is natural to think of a "bosonic counterpart" of the exclusion 
process. In fact, such a process was introduced in [H] and [TO] as the dual of the 
Brownian momentum process, a stochastic model of heat conduction (similar models 
of heat conduction were introduced in [I] and [8], see also [5] for the study of the 
structure function in a natural asymmetric version). 

In the present paper we analyze this "bosonic counterpart" of the exclusion pro- 
cess. We will call this process (as will be motivated by a Poisson clock representation) 
the "symmetric inclusion process" (SIP). In the SIP, jumps are performed according 
to independent random walks, and on top of that particles "invite" other particles to 
join their site (inclusion). For this process we prove the analogue of the comparison 
inequality for the symmetric exclusion process. From the comparison inequality, using 
the knowledge of the stationary measure and the self-duality property of the process, 
we deduce a series of correlation inequalities. Again, in going from exclusion to inclu- 
sion process the correlations turn from negative to positive. We remark however that 
these positive correlation inequalities are different from the ordinary preservation of 
positive correlations for monotone processes [12], because the SIP is not a monotone 
process. Since the SIP is dual to the heat conduction model it is immediate to extend 
those correlation inequalities to the Brownian momentum process and the Brownian 
energy process. 

We also introduce the non-equilibrium versions of the SIP, i.e., we consider the 
boundary driven version of SIP. In this case, for appropriate choice of the boundary 
generators, we prove duality of the process to a SIP model with absorbing boundary 
condition. We then deduce a correlation inequality, explaining and generalizing the 
positivity of the covariance in the non-equilibrium steady state of the heat conduction 
model in [9]. 

All the results will be stated in the context of a family of SIP(m) models, which are 
labeled by parameter m e N. As the SEP model can be generalized to the situation 
where there are at most n G N particles per site (this corresponds to a quantum spin 
chain with SU(2) symmetry and spin value j = n/2), in the same way the SIP model 
can be extended to represent the situation of a quantum spin chain with SU(1,1) 
symmetry and spin value k = m/4 [9]. 

The paper is organized as follows. In Section 2 we define the SIP(m) process, re- 
stricting to a context where its existence can be immediately established. The main 
comparison inequality, which allows to compare SIP walkers to independent walkers 
(by a suitable generalization of Liggett comparison inequality) is proved in Section 3. 
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Correlation inequalities for the SIP(m) process that can be deduced from the com- 
parison inequality are proved in Section 5 (the necessary knowledge of the stationary 
measure and the self-duality property are presented in Section 4). In particular, in 
Section 5 it is proved that when the SlP(ra) process is started from its stationary 
measure then correlations are always positive, while when the process is initialized 
with a general product measure then positivity of correlations is recovered in the long 
time limit. Further correlation inequalities for systems similar to the SlP(ra) process 
are discussed in the subsequent Sections. Attractive interaction (the SEP(n), which 
generalize the standard SEP) is presented in Section 6. Some interacting diffusions 
dual to the SIP(m) process are studied in Section 7. Finally the boundary driven 
SIP(m) process is analyzed in Section 8. 

2 Definition 

In the whole of the paper, S will denote or a finite set, or S = Z d . Next, p(x, y) 
denotes an irreducible (discrete-time) symmetric random walk transition probability 
on S, i.e., p(x, y) = p(y, x) > 0, p(x, y) = 1, and p(x, x) = 0. In the case S = Z d , 
we suppose furthermore that p(x,y) is finite range and translation invariant, i.e., 
p(x, y) = n(y — x), and there exists R > such that p(x, y) = for \x — y\ > R. This 
assumption for the infinite-volume case avoids technical problems for the existence 
of the SIP(m) which for the subject of this paper are irrelevant. The proof of 
existence of the SIP(m) in our infinite- volume context (with the process started 
from a "tempered" initial configuration, i.e. rj(y) < \\y\\ k for some k and for all y) 
follows from self-duality, along the lines of [6], Chapter 2. 

The symmetric inclusion process with parameter m £ (0, oo) associated to the 
transition kernel p is the Markov process on O := N s with generator defined on the 
core of local functions by 

Lf(v)= ^p(x, 2 /)2r /x (m + 2^)(/(r / ^)-/(r / )) (2.1) 

x,y£S 

where, for r\ £ Q, t] x,y denotes the configuration obtained from rj by removing one 
particle from x and putting it at y. 

In [S], for m = 1 this model was introduced as the dual of a model of heat 
conduction, the so-called Brownian momentum process, see also [TU], and [1] for 
generalized and or similar models of heat conduction. 

The process with generator (12. ip can be interpreted as follows. Every particle has 
two exponential clocks: one clock -the so-called random walk clock- has rate 2m, the 
other clock -the so-called inclusion clock- has rate 4. When the random walk clock 
of a particle at site x £ S rings, the particle performs a random walk jump with 
probability p(x, y) to site y £ S. When the inclusion process clock rings at site y £ S, 
with probability p(y, x) = p(x, y) a particle from site x £ S is selected and joins site 

V- 
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From this interpretation, we see that besides jumps of a system of independent 
random walkers, this system of particles has the tendency to bring particles together at 
the same site (inclusion), and can therefore be thought of as a "bosonic" counterpart 
of the symmetric exclusion process. 

To make the analogy with the exclusion process even more transparent, in an 
exclusion process with at most n particles (n G N) per site (notation SEP(n)), the 
jump rate is r]i(n — rjj)p(i,j). Apart from a global factor 4, the SIP(m) is obtained 
by changing the minus into a plus and choosing n = m/2. 

Notice that the rates in (12.1 ft are increasing both in the number of particles of the 
departure and in the number of particles of the arrival site (the rate is p(x, y)2r) x (m + 
2r) y ) for a particle to jump from x to y). Therefore, by the necessary and sufficient 
conditions of [IT], Theorem 2.21, the SIP is not a monotone process. It is also easy 
to see that due to the attraction between particles in the SIP, there cannot be a 
coupling that preserves the order of configurations, i.e., in any coupling starting from 
an unequal ordered pair of configurations, the order will be lost in the course of time 
with positive probability. 



2.1 Assumptions on the transition probability kernel 

In this section we introduce the assumptions that we need to prove the positivity 
of correlations of stationary measures obtained as limits of general initial product 
measures (see later for precise definitions). This assumptions are only relevant in the 
infinite volume case S = Z d and they are indeed satisfied in the context of finite-range 
translation-invariant underlying random walk kernel p(x, y) = ir(y — x). However, all 
our results on correlation inequalities for stationary measures depend only on one or 
both of the assumptions below, i.e., if on more general graphs, or on Z d with more 
general p(x,y), existence of SIP(m) would be established, then the corresponding 
correlation inequalities hold under one or both of the assumptions Al, A2 below. 

We define the associated continuous-time random walk transition probabilities of 
random walk jumping at rate 2m: 

Mx,y) = f2^^e- 2m V n) (x,y) (2.2) 

n=0 

where p^ denotes the n th power of the transition matrix p. Denote by the 
probability measure on path space associated to two independent random walkers 
X t , Y t started at x, y and jumping according to ( 12. 2p and by Pf,y P ^ the corresponding 
probability for two SIP walkers X t , Y t jumping with the rates of generator (12. ip . 

We consider two assumptions 
- Assumption (Al) 
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Assumption (A2) 

lim sup Pff/ (m) (X' t = y t ') = (2.4) 



The assumption (Al) amounts to requiring that for large t > 0, two independent 
random walkers walking according to the continuous time random walk probability 
(12. 2p will be at the same place with vanishing probability. The assumption (Al) 
follows immediately if we have 

lim supp t (x, y) = (2.5) 

x,y 

since then 

lim sup P' i ? y(m) (X t = Y t ) = lim sup p t (x, u)p t {y, u) = lim supp 2t (x,y) = 

x,V x,y x,y 

(2.6) 

Notice also that, by simple rescaling of time, (Al) holds for all m > as soon as it 
holds for some m > 0. 

Assumption (A2) guarantees that two walkers evolving with the SIP dynamic will 
be typically at different positions at large times. Notice that in the case we consider, 
i.e., the translation invariant finite-range case S = Z d , p(x, y) = p(0, y—x) =: ir(y—x), 
this is automatically satisfied, as the difference walk X t — Y t of two SIP particles is 
a random walk Z t on Z d with generator 

L z f(z) = 8tt(*)(/(0) - /(*)) + ^4mvr(y)(/(z + y) - f(z)) (2.7) 

y 

which is clearly not positive recurrent. 

Assumption (A2) implies that any finite number of SIP particles will eventually 
be at different locations. This is made precise in Lemma [TJ in section 



3 Comparison of the SIP with independent ran- 
dom walks 

We will first consider the SIP process with a finite number of particle in subsection 
13.11 and then state the comparison inequality in subsection 13.21 

3.1 The finite SIP 

If we start the SIP with n particles at positions xi, . . . ,x n G S, we can keep track 
of the labels of the particles. This gives then a continuous-time Markov chain on S n 
with generator 

n / n \ 

c n f( Xl ,...,x n ) = J2T,M^y)[rri + 2j2Hy-^))(f^ Xi ' y )-f(^) 

i=l y€S \ j=l J 

= £ ltn f(x)+£ 2>n f(x) (3.1) 
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where x XuV denotes the n-tuple (xi, . . . , y, x i+ i, . . . , x n ). Further, £i )n , resp. £ 2 ,n 
denote the random walk resp. inclusion part of the generator and are defined as follows 

n 

c 1 , n f(x 1 , ..., Xn ) = J2Yl 2m P^ y)(f( xX " y ) - /(*)) ( 3 - 2 ) 

i=l yes 

n n 

C 2 , n f{ Xl , ...,*„) = Y,Y d Ap{x ii x j ){f{x B ^) - f{x)) (3.3) 

i=l j=l 

3.2 Comparison inequality 

From the description above, it is intuitively clear that in the SIP, particle tend to be 
less spread out than in a system of independent random walkers. Theorem [T] below 
formalizes this intuition and is the analogue of a comparison inequality of the SEP 
([H], Chapter VIII, Proposition 1.7). 

To formulate it, we need the notion of a positive definite function. A function 
/ : S x S — > R is called positive definite if for all /3 : S — > R such that Y2 X < 00 

X)/(x,j/)/3(x)/3(y)>0 

x,y 

A function / : S n — > R is called positive definite if it is positive definite in every pair 
of variables. 

We first introduce a slightly more general generator with parameters a > 0, b G R 
that includes both process of exclusion and inclusion type. 

n / n \ 

C?f(x u ...,x n ) = £ X>(x„ y)la + bJ2Hy = *i) (/0=*" w ) - /(*)) (3-4) 
i=l yes \ j=i J 

SO 

pa, b pa I pb 

where 

n 

ClJXx u ...,x n ) = a ^ ZX^ //)(•/'(•''"•") " /(*)) ( 3 - 5 ) 

i=l yeS 

is the independent random walk part (random walks jumping at rate a) and 

n n 

) )(/(.''•'•'••'•' ) - /(.,•)) (3.6) 

i=l i=i 

is the "clumping" part, i.e., when b < clumping is discouraged, and b > clumping 
is favored. 

We call T£ ,b (t) the semigroup on functions / : S n — > R associated to the generator 
(13. 4p . and £/^(£) the semigroup of a system of independent continuous-time random 
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walkers (jumping at rate a), i.e., the semigroup associated to the generator L\ n in 
(13. 5p . Notice that when b < 0, T£' b (t) is not always a Markov semigroup. However, 
for the applications of negative b, we have in mind generalized exclusion process 
(see Section 6) in which case a/b is an integer and in this case T£ ,b (t) is a Markov 
semigroup. 

Theorem 1. Let f : S n — > R be positive definite and symmetric. Then we have for 
b>0 

K(t)f < T^(t)f (3.8) 
and for b < 0, if (T a,b (t)) t >o is a Markov semigroup, we have 

K(t)f > T^ b {t)f (3.9) 

Proof. The proof follows the proof in [2], but for the sake of self-constistency we 
prefer to give it explicitely. Suppose b > 0. 

Start with the decomposition (13. ip and use the symmetry of p(x, y) and / to write 
(C a /f-ClJ)(x) = (C b f)(x) 



bY,Y,l' U '- x ^f( :r '' '''I'- /(*)) 



i=i j=i 

n n 



/.(./', .o)(/(.r'""- ) + /(,••'-' ) - 2f(x)) 



2 

i=i j=i 
b n n 

i=l j=l 

X ^ ^ f(.%ll • • • ) 1) • • • ) «Cj — 1 j S'j + l, • - - , 3?n) i$Xi,U 0~Xj,u) i$Xi,V 0~Xj,v) 

u.v 

> (3.10) 

where in the last step we used that / is positive definite. 

Since U°(t) is the semigroup of independent walks, it maps positive definite func- 
tions into positive definite functions, and so we have 

{c n u a n {t)f - c a lt Mt)f) = cl n u a n (t)f > 

We can then use the variation of constants formula 

T:> b (t)f - U a n (t)f = fds T?{t - s) {Cl n U a n {s)f) > (3.11) 
Jo 

and remember that T° ,b (t) is a Markov semigroup which therefore maps non-negative 
functions into non-negative functions. 

The proof for b < 0, under the assumption that T^ b {t) is a Markov semigroup is 
identical. □ 
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4 Stationary measures and self-duality for the SIP(m) 

The stationary measures of SIP(m) are product measures of "discrete gamma distri- 
butions" 

where for n > 

with < A < 1 a parameter, T(r) the gamma-function and 



; A,m 



x m/2 



A 



Notice that for m = 2, ^™ is a geometric distribution (starting from zero), i.e., 
z/^(n) =A"(l-A),nGN and for m/2 an integer v™ is negative binomial distribution 
NB(m/2, A). Moreover, the measures v m have the following convolution property 

A A = ^A I 4 - 2 ) 

where * denotes convolution, i.e., a sample from i/™ * v\ is obtained by site- wise 
addition of a sample from v™ and an independent sample from v\. 

The SIP(m) process is self-dual [10] with duality functions given by D(£,rj) = 
X[ x d{£,x,r)x), with 

7| T (-) 

d(k, I) = - — . u ' . (4.3) 
1 ; (I- k)\T(f + k) K J 

where k < I. Self-duality means that 

K P{m) m Vt) = Ef P{m) D(C h rj) (4.4) 

where E^ /P ^ denotes expectation in the SIP process started from the configuration 
V- 

The relation between the polynomials D and the measure v™ reads 

D(t,r,Wfa)=(j^y (4.5) 

as follows from a simple computation using the definition of the T-function, T(r) = 
J °° x r ~ 1 e~ x dx. 

From conservation of particles in the dual process, we see that self-duality and 
the relation ( 14. 5 p gives stationarity of the measure v\. 

The relation ( 14. 5 p can be generalized to "local stationary measure", i.e. the prod- 
uct measures that are obtained from the stationary measure (14. ip by allowing a site- 
dependent parameter. More precisely, given 

A: S-> [0,1) 
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we define the local stationary measure associated to the profile A by 

^=®xesv% B) {dr h ) (4.6) 

For x±,...,x n G S we denote by Yli=i ^ the particle configuration (eN s obtained 
by putting a particles at locations x i} i.e., £(x) = Y^i=i H x i = x )- We then have the 
following relation between the duality functions and the local stationary measures 

„ / n \ n 

J D [Y^S^A n(d V ) = l[p(x t ) (4.7) 



where 



\(Xi 



1 - X(xi) 



For a constant profile X(x) = A ,Vx G S, we recover (14. 5p . 

By Lemma [T] below, in the case S = 1* d and translation invariant finite-range 
p(x,y), any number of dual particles in the SIP{m) will eventually diffuse away to 
infinity. From that it is easy to deduce that the measures v\ are extremal invariant. 
To see this, we denote for two finite particle configurations £ _L £', if their supports are 
disjoint, i.e., there are no site x G S where there are £ and £' particles. If £ _L £' then 
D(£ + £',??) = -D(£, ? ?)-D(£ / , v)- Since at large t > 0, assumption (A2) implies that, 
in the SIP started with a finite number of particles, particles are with probability 
close to one at different locations (see Lemma [1] for a proof of this), we have that 
for £' a fixed configuration, the event £ t _L £' has probability close to one as t — > oo. 
Therefore 

hm / E s v Ip ^(D(^ Vt ))D(e,V>x(d V ) 



hmEf P(m) / D(C t , V )D(e,V>x(dv) 

t— >oo * J 

hmEf P(m) / D(£ t ^)D(£^)/(6 JL C>x(dv) 
hm 

t— >oo 



which shows that time-dependent correlations of (linear combinations of) -D(£, •) poly- 
nomials decay in the course of time to zero, and hence, by standard arguments, v\ is 
mixing and thus ergodic. 



9 



5 Correlation inequalities in the SIP(m) 



For a probability measure \i on the configuration space N s , we denote its "duality 
moment function" : S n — > K by 



K^xi, . . . , x n ) = J D ^2 tisuvj K d v) (5-i; 



If fj, = ur is a local stationary measure with profile A, then 



n 



K^(xi, . . . ,x n ) = Y\_p(xi) (5.2) 



i=l 



which is clearly positive definite and symmetric. We can therefore apply Theorem [T] 
and obtain the following result. 



Proposition 1. For all t > 0, for all profiles A : S — > [0, 1) and for all xi, . . . , x n G S 

we have 



K^st (xi, ■ ■ ■ , x n ) > JJ (x^ (5.4) 



i=l 



where S t denotes the semigroup of the SIP(m) process. In particular, when the 
SIP{m) is started from u^, the random variables {r)t(x),x G S} are positively corre- 
lated, i.e., for (x,y) G S x S 



^ IP{m) (Vt(x) Vt (y))i^(d v ) > / E s v Ip ^(r) t (x))^(dr]) / Ef p W (r} t (y)) u^dr,) 



Proof. Denote by Ef^f™] expectation in the SIP(m) process started with n particles 
at positions (x±, . . . ,x n ), by ^ IRW ( m ) expectation in the process of independent ran- 
dom walkers (jumping at rate 2m) and ~E RW ( m ) a single random walker expectation. 
We then have the following chain of inequalities, which is obtained by using sequen- 
tially the following: self-duality property (I4.4p . the comparison inequality (13. 8p . the 
relation between the measure and the duality function D (14. 7p . the independence 
between random walkers, the fact that a single SIP particle moves as a continuous 
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time random walk, and finally again self-duality ( 14. 4|) 

E siP {m)D (|>^ u^dn) 

(n N 

n 

i=i 

n „ 

= II / C M (D(S Xi{t) , V )) u^dr,) 



i=l 



n „ 

= n K p{m) ( D (^))v- x (d V ) (5.5) 
i=i j 

□ 

This proposition shows that starting from a local stationary measure v^, the den- 
sity profile pt(x) = E, x W ^ p t (x) predicts (by duality) correctly the density at time 
t > but the true measure at time t > 0, v^St , lies above (in the sense of expectations 
of .D-functions) the product measure with density profile Pt{x). 

From the analogy with the SEP emphasized above, one could think that (15.41) 
extends to the case when the SIP process is started from a general product measure. 
However, for general probability measures /i on Q, the duality moment function : 
S n — > R defined in (15 .ip is not necessarily positive definite (as is the case for the 
special product measures z^-), since we do not have the equality D (X)£=i S Xi , r/) = 
Q™ =1 D(5 Xi ,r]) in general. Notice that this problem does not appear in the context of 
the standard SEP, as for that model, the self-duality functions are 

(n \ n n 

j2 v ) = n ^ = n DsEp i 5 ^) 
i=l ) i=l i=l 

and hence automatically, for any measure p, the function is positive definite in 
that model. 

If however all Xi are different, we have £)(X^Li $xi,v) = YYi=i D(5 Xt ,ri). For every 
probability measure p on Q, the function ^ : S n — > 1R defined by 

/n 
Y[D(5 XiJ 7])fi(dr]) (5.6) 
i=i 
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is clearly positive definite. This, together with the fact that under assumption (A2), 
a finite number of SIP(m) particles diffuse and therefore eventually will be typically 
at different positions, suggests that in a stationary measure, the variables r\ Xi are 
positively correlated. 

To state this result we introduce the class of probability measures with uniform 
finite moments 

V f =: {/j, : Vn £ N, sup / D(£, rj)K d v) =■ M? < 00} (5.7) 

\t\=nJ 

For a sequence of measures fi n £ Vf, and fj, £ Vf, we define that fi n — > fi if for all £ 
finite particle configuration, 



lim / D(£,r])v n (dr}) = / D(£,r])[i(dr]) 



We can then formulate our next result. 

Proposition 2. Assume (Al) and (A2). Let v £ Vf be a product measure. Let S(t) 
denote the semigroup of the SIP(m). Suppose that 

\x = lim vS{t n ) (5.9) 

n— ¥00 

for a subsequence t n | 00. Then we have \i £ Vf, fi is invariant and 

n 

K„(x 1 ,...,x n )>Y[K l ,(x i ) (5.10) 

i=i 

Proof. First, by duality we have, referring to the definition of Vf, for all t > 0, 



which shows that both vS{t n ) and \i are elements of Vf. The invariance of \i follows 
from duality, v £ Vf and Lemma 1.26 in [H], chapter V. 

To proceed with the proof of the proposition, we start with the following lemma, 
which ensures that, under condition (A2), any number of SIP(m) particles will even- 
tually be at different locations. 

Lemma 1. Assume (A2). Start the finite SIP(m) with particles at locations {x\, . . . , x n }, 
then 

Mm P£22 (3i ^ J : Mt) = X i(t)) = (5.12) 
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Proof. We give the proof for m = 1. The general case is a straightforward extension. 
Put 77 := Ym=i $xi- Using self-duality we can write 

Pf P(1) (3^j:X,(t)=X,(t)) < J2 F v IPil) (Vt^)-Vt(z)>l) 

z 

< E E f p(1) (*)-^)) 



±J2K m (D(26 z , Vt )) 

z 
z 



< 



= 3 E E E f/ (1) = **) w = 

< 3n 2 supFfi P ^(X t = y t ) 



(5.13) 



where in the last step we used the symmetry of the transition probabilities of the 
SIP(1) (with two particles). □ 

We now proceed with the proof of the proposition. For xi, . . . , x n G S we define 



1=1 



1=1 



A(xi, . . .,X n ,7]) 



(5.14) 



We have that A(x±, . . . , x n , 77) = if all Xj are different, i.e., if \{xi, . . . , x n }\ = n. 
Since by assumption (A2) and Lemma [TJ the probability that two SIP(m) walkers 
out of a finite number n of them occupy the same position, i.e. X{(t) = Xj(t) for some 
i 7^ j, vanishes in the limit t -> 00, we conclude, using v G Vf, for any Xi, . . . , x n G S, 

lim / Ef i / i P SA(X 1 (t),...,X n (t),r / )z/(^) = (5.15) 

Moreover from the comparison inequality (13. 8p we have, using the notation (15. 6p 
Kl^^X^tl . . . , X n (t)) > EjfWW), . . . , X n (t)) 

i=l 

n „ 

= n E r (m) / ^fe( t) ^)z/(^)+6(t) (5.16) 



i=l 



where e(t) — > as t — > 00 by assumption (Al), i.e., for large t > 0, independent 
random walkers are at different locations with probability close to one. Therefore, 



13 



using the definition (15 .9|) . the self-duality property (I4.4p . the equation (15.151) . the 
equation (15.161) . and taking limits along the subsequence t n we have 



Ku,{pC\ i • • • i •En) 



i=i 



C) 



77 ) ^(dr/) 



= limE^M^x^),...,^)) 

t— >oo ' ' 

n „ 

> t hmnEr M / D(5 Xi(t)l r ] )u(d V ) 
i=i 



l[Kp(xi) (5.17) 

□ 



i=l 



6 Correlation inequalities in the SEP(n) 

We now consider the application of the generalized Liggett inequality for negative b. 
The SEP(n) is the Markov process on Q = {0, 1, . . . , n} s with generator 

Lf(v)= ^v(x)(n-T,(y))p(x,y)(f(Tf*)-f(T,)) (6.1) 

x,y€S 

The stationary measures of this process are products of binomial distributions, i.e., 
for pe [0,1], 

v p = ®xesBin(n, p) (6.2) 

Similar to the case of the inclusion process, for a profile p : S — > [0, 1] we define the 
local stationary measure 

Vp = <3 x &sBin(n,p(x)) 
The duality functions for self-duality are given by (see [TO] ) 



D{^v) = Hd(^rj x ) (6.3) 

X 

for £ G Q a configuration with finitely many particles (at most n per site) and with 

d(M) = 4l (M 

The relation between the duality functions and the local stationary measures is, as 
usual, i.e., for £ = XX=i ^< e ^ (*- e -' a ^ mos t n particles per site), and p a profile: 

„ n 

J D(^r ] )^(d V ) = l[p(x l ) (6.5) 
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We define, for a probability measure fi on f2, its duality moment function 

fj,(drj) (6.6) 



K^xx, . . . , x n ) = J D ^2 



The following proposition is then the analogue of Proposition [T] in this context 
(with inequality in the other direction since b < 0). 

Proposition 3. For J) : S — > [0, 1] a density profile and t > 0, 

n 

K^ St {xi, ...,x n )< Y[K^ St (xi) (6.8) 

i=i 

In particular, for starting from Vp, the variables {r]t{x) are negatively corre- 

lated. 

7 Correlation inequalities for some interacting dif- 
fusions 

7.1 The Brownian Momentum Process 

The Brownian momentum process is a system of interacting diffusions, initially in- 
troduced as a model of heat conduction in and analyzed via duality in [SJ. It is 
defined as a Markov process on X = ~R S via the formal generator on local functions: 



L B Mpfiv) = ( p( x >y) f 

\x,yeS ^ 



d d N 



drj y 8% 



-V5- /to) 



The variable r] x has to be thought of as momentum of an "oscillator" associated to 
the site x G S. The local kinetic energy r\ 2 x has to be thought of as the analogue of 
the number of particles at site x in the SIP(m) with m — 1. The expectation of t] x 
is interpreted as the local temperature at x. 
Defining the polynomials 



D(n,z) 



z 2n 



(2n- 1)!! 



we have the duality function D(£, •) defined on X and indexed by finite particle 
configurations £ G N 5 , J2 x £ x < oo: 

D(t,7i) = Y[D(Z x ,r, x ) (7.2) 

In [pj, [TO], we proved the duality relation 

E* MP (£>(£, Vt)) = Ef P(1) (£>(& r/)) (7.3) 



15 



As before, for x\, . . . , x n G S we denote by Y^=i t ne particle configuration obtained 
by putting a particle at each Xj. 

Let fi be a product of Gaussian measures on X, with site-dependent variance, i.e., 
for a function p : S — > [0, oo), we define 

Hp = ®xesv p (x)(drix) (7.4) 

where 

e -»?i/2p(s) 

v p {x){dr) x ) = , dr} x 
v /2vrp(x) 

is the Gaussian measure on M with mean zero and variance p(x). Then we have 

{Xi) (7.5) 



// n \ n 

\i=i j i=i 



From this expression, it is obvious that the map 

S n ^R:(x 1 ,...,x n )^ [ D I > ]6 Xi ,T}\n P (dr}) (7.6) 



is positive definite. Therefore, combining the duality property between BMP process 
and SIP(m) process, (17.31) . with Theorem [1] we have the inequality 



J E* MP D (X>^*J M*?) 



KZS (j[jD(S Xiit) ,v)^(dv)] 

n 

ri E r (m) p(^(t)) 

8=1 

n „ 

n/Efr« (2?fe (t) ^))A* P W 

n „ 

n / (D (5 Xi} r)t)) ^ P {dr}) (7.7) 

i=i J 
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which is the analogue of Proposition [JJ for the BMP process. 

In words, it means that the "non-equilibrium temperature profile" is above the 
temperature profile predicted from the discrete diffusion equation. It also implies 
that the variables {r)% :i£5} are positively correlated under the measure (p p )t for 
all choices of p, t > 0. 

More precisely, if we denote 

Pt (x) = E* w ™p(X t ) 

then we have that rf x at time t has expectation pt{x) when the starting measure is p p 
(since a single particle in the SIP(1) moves as a continuous time random walk). The 
correlation inequality for the BMP which we just derived shows that the true measure 
at time t > when started from a product of Gaussian measures lies stochastically 
above the Gaussian product measure with mean zero and variance pt{x). 

Similarly, we obtain an analogous correlation inequality for the BMP for a measure 
obtained cis du limit of product measures. We define 

V f (X) = {p : Vn e N : sup / D(£,r})p(dr)) < 00} 

\(\=nJ 

Proposition 4. Assume (Al) and (A2). Suppose v e Vf(X) is a product measure 
and p is a limit point of the set {vS(t) : t > 0}, where S(t) denotes the semigroup of 
the BMP process. Then we have the inequality 

n 

Kn(xi, ...,x n )> Y\_K^(xi) 

*=i 

7.2 The Brownian Energy Process 

The Brownian energy process with parameter m > (notation BEP(m)) is intro- 
duced in [10] as the process on state space X = [0, 00) s , with generator 

L ( 7 - 9 ) 

x,yeS 

with 

W = *»* (£ - - 2K* - %) - £) /W (7.10) 

This process is dual to the SIP(m) in the following sense. Define, for £ e N 5 a finite 
particle configuration, and r] E X the polynomials 

xg5 
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with, for k e N, y e [0, oo) 



r (—) 

d(k, y) = y k , M 2 \ (7-12) 



then we have 



E BEP(m) D ^ Vt) = E 8IPim) D fa v) (7 . 13) 

As a consequence, extremal invariant measure of the BEP(m) are products of re- 
distributions with shape parameters m/2 and scale parameter 9 > 0: 

u d (drj) = <S> xe s^e(dT]x) (7.14) 

with 

u e (dz) = . 1 , r zT-i e ~Vg (7.15) 

n ; ^/2r(s) v ; 

Similarly we define the local stationary measures 

*fr= ®xeSVd( x) (dr]x) (7.16) 

with : S 1 — >■ [0, oo), and the duality moment function of a probability measure fi on 
X: 

K ll (x 1 ,...,x n )= f D ^6^,7) ]fi(dri) (7.17) 



As a consequence of the correlation inequalities derived for the SIP(m), we derive 
the following. 

Proposition 5. i. For all 9 : S 1 — >■ [0, oo), t > 0, and xi, . . . ,x n G S we have 

n 

K^BEPim) (X 1 , . , X n ) > Yl K^ s BBP(m)fl(Xi) (7.19) 

i=l 

2. If for some product measure v on X with finite moments, and a sequence of 
t n 'f oo the limit 

/i= lim vS BEP{rn \t n ) 

n—too 

exists, then 

n 

K fl {x 1 ,...,x n )>Y[K fl {x i ) (7.20) 

i=i 



8 The boundary driven SIP(m) 

In this section we consider the non-equilibrium one-dimensional model that is ob- 
tained by considering particle reservoirs attached to the first and last sites of the 
chain. We will show that, if one requires reversibility w.r.t. the measure v™ and 
duality with absorbing boundaries, this uniquely fixes the birth and death rates at 
the boundaries. 
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8.1 Duality for the the boundary driven SIP(m) 



The generator of the boundary driven SIP(m) on a chain {1, . . . , N} driven at the 
end points, reads 

C = C\ + Cn + C-buik (8.1) 

where Cbuik denotes the SIP(m) generator, with nearest neighbor random walk as 
underlying kernel, i.e., 

Cbuikfiv) = E 2^(m+2 %+1 ) {f(v x ' x+1 ) ~ f(v))+2v x+ i(m+2 Vx ) (/(rf +1 >*) - f(rj)) 
xe{i,...,N~i} 

(8.2) 

and where C\,Cn are birth and death processes on the first and iV-th variable re- 
spectively, i.e., 

Cif(v) = d L (vi)(f(v ~ *i) - f(v)) + b L (vi)(f(v + *o - m) 

and 

C N f(r]) = d R (rj N )(f(r] - 5 N ) - f( V )) + b R ( VN )(f( V + 5 N ) - f( V )) 

These generators model contact with respectively the left and right particle reservoir. 

The rates d,L, b^, d R , b R are chosen such that detailed balance is satisfied w.r.t. the 
measure v™, with A = X L for d>L,b L , and X = \r for d R ,b R . More precisely, this 
means that these rates satisfy 

b a (k)uZ(k) = d a (k + l)vZ{k + 1) (8.3) 

for a G {L, R}. 

To state our duality result, we consider functions T>(£,rj) indexed by particle 
configurations £ on {0, . . . , N + 1} defined by 

V(^v) = P^D(Z {1 ,..., Nh r ] )p^ (8.4) 
where p a = p\ a = \ a /{\ — A Q ), and where we remember that 

n! T (f ) 



D{k,n) 



(n-fc)lr(f + A;) 



is the duality function for the SIP(m). I.e., for the "normal" sites {1,...,N} we 
simply have the old duality functions, and for the "extra added" sites {0, iV + 1} we 
have the expectation of the duality function over the measure v™. 

We now want duality to hold with duality functions £>, and with a dual process 
that behaves in the bulk as the SIP(m), and which has absorbing boundaries at 
{0, N + 1}. More precisely, we want the generator of the dual process to be 

£ = C-bulk + £>i + C-n (8.5) 
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with C bulk given by (JS2J, and 

^/(o = ^ - /(e)) 

for £ G N't ' 1 ---' JV+1 }. The duality relation then reads, as usual, 

(CV^.))( V )=(CV(., V ))(0 (8.6) 
Since self-duality is satisfied for the bulk generator with the choice (18. 4p . i.e., since 

(18. 6 p will be satisfied if we have the following relations at the boundaries: for all 
k < n: 

b a (n)(D(k, n+1)- D(k, n)) + d a (n)(D(k, n - 1) - D(k, ri)) 
= k{D(k-l,n)p a - D(k,n)) (8.7) 

where a G {L, R}. 

From detailed balance (18. 3p we obtain 

^^KfT^i)^"- 1 ) (8 ' 8 > 

Working out (18. 7p gives, using ( 14. 3p . 



( "t 1 , -l)+da(n) (—-1 



n+1 — k I \ n 



which simplifies to 



b a (n) _ dM = I (f + k-l)p Q _ ^ 
n + 1 — k n \ n — k+ 1 



(8.10) 



Choosing 

and by the detailed balance condition (18. 8 p 



TL 

da(n) = - r- (8.1i; 

1 — A Q 



W">= Gr + ")i^r < 8 - 12 > 
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it is then an easy computation to see that (18.71) is satisfied with the choices (18. lip . 
(I8.12p . Indeed, (I8.10p reduces to the simple identity 

fm \ / A \ 1 1_ f + k- 1 / A \ 

V 2 +n )\l-Xjn+l-k~l-X~ n + l- k \l- XJ ~ 

We remark that the requirement of detailed balance alone is not sufficient to fix the 
rates uniquely. However, the additional duality constraint ( 18. 7p does fix the rates to 
the unique expression given by (18. lip and (I8.12p . 

As a consequence of duality with duality functions (18. 4p . we have that the bound- 
ary driven SIP(m) with generator ( 18. ip has a unique stationary measure ^l,r for 
which expectations of the polynomials D(^,n) are given in terms of absorption prob- 
abilities: 

J t— ¥00 

t— >oo 

= E PlPr^Z (Zoo = k6 + l6 N+1 ) (8.13) 

k,l:k+l=\£\ 

Here, denotes expectation in the dual process (which is absorbing at {0, N + 1}) 
starting from £. In particular, since a single SIP(m) particle performs continuous 
time simple random walk (at rate 2m) we have a linear density profile, i.e., 

J D(5 t ,n)^ R (dn) = p L (l- + Prj^ (8.14) 



8.2 Correlation inequality for the boundary driven SIP(m) 

For xi,...,x n G {1, . . . , N} let us denote by (Xi(t), . . . , X n {t)) the positions of par- 
ticles at time t evolving according to the SIP(m) with absorbing boundary sites at 
{0, iV + 1}, i.e., according to the generator (18. 5p . and initially at positions x\, . . . , x n . 
Let (Yi(t), . . . , Y n (t)) denote the positions at time t of independent random walk- 
ers (jumping at rate 2m) absorbed (at rate 1) at {0, N + 1}, initially at positions 
X\,...,x n . Since the absorption parts of the generators of (Xi(t), . . . ,X n (t)) and 
(Yi(t), . . . , Y n (t)) are the same, we have the same inequality for expectations of pos- 
itive definite functions as in Theorem [TJ Therefore, we have the following result on 
positivity of correlations in the stationary state. This has once more to be com- 
pared to the analogous situation of the boundary driven exclusion process, where the 
stationary covariances of site-occupations are negative. 

Proposition 6. Let Hl,r denote the unique stationary measure of the process with 
generator ( 18. ip . Let x%, . . . ,x n G {1, . . . , N}, then we have 

J D (X>^1 »L, R (dri) >f[J D(6 Xi ,r))v L>R (dr)) (8.16) 
In particular, n x ,x G {1, . . . , N} are positively correlated under the measure ^l,r- 
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Proof. Start from the measure v™. Define the map {0, . . . , N + l} n — > R: 

// n \ n 

V i=i J i=i 

where p(x) = for x 6 {1, ... , iV} and p(0) = /?£,, p(iV + 1) = p^. This is clearly 
positive definite. Therefore, for Xi, . . . , x n G {1, . . . , N}, we have 

= lim J Ef^M-^ f P^fr^.ij) J v?(dr)) 



n 

n „ 

= II / DiS^^^Ridrj) (8.18) 

i=l ^ 

where we denoted E s/p ( m )> a6s for expectation over SIP(m) particles absorbed at 
{0,N + 1}, and E IRW( - m ^ abs for expectation over a system of independent random 
walkers (jumping at rate 2m) absorbed (at rate 1) at {0, N + 1}. □ 

Remark 1. 1. Proposition® is in agreement with the findings of JPJ/, where the 
covariance ofr]i,r]j in the measure Pl,r was computed explicitly, and turned out 
to be positive. 

2. For the nearest neighbor SEP on {1, . . . , N} driven at the boundaries, we have 
self-duality with absorption of dual particles at {0, N + 1} and duality function 

(n \ n 

i=l / t=l 

where r] := Pl,Vn+i = Pr- Since for SEP particles we have the comparison 
inequality of Liggett, we have as an analogue of ( I8.16P in the SEP context, 

/n n „ 

i=l i=l J 

i.e., 7] Xi are negatively correlated. The same holds for the non- equilibrium 
SEP{n) driven by appropriate boundary generators. This is in agreement with 
the results in JW$ . where the two-point function of the measure pl,r is com- 
puted, and with the work of ^j, where some multiple correlations are explicitly 
computed. 
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3. We expect the KMP-model, a model of heat conduction introduced and studied in 
ITSf to also have positive correlations. Indeed, the KMP and the BEP(2) model 
are related by a so-called instantaneous thermalization limit 110]. Therefore, it 
is natural to think that similar correlation inequalities should hold for the KMP 
as we have derived for the BEP. The limit to obtain the KMP from the BEP 
is however difficult to perform on the level of the n-particle representation and 
it is thus not clear (to us) how to prove that the KMP preserves the positive 
correlation structure of the BEP. A positive hint in this direction comes from 
the explicit expression of the two point function which has been computed for 
the KMP in the non- equilibrium context in J$j. 
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